// Volumatric flux
surfaceScalarField& phi = fluid.phi();
